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Magnetic reconnection is a ubiquitous plasma process that converts magnetic energy 
into thermal and kinetic energy. Of particular interest is the production of non-thermal 
particles, in which a fraction of the plasma population is driven to energies much larger 
than that found in the ambient medium. Reconnection is thought to be an important driver 
of such particles in phenomena such as gamma ray bursts 1,2 , stellar and solar flares 3 , and 
magnetospheric storms 4 . Recent observations of solar flares reveal the remarkable efficiency 
of electron acceleration: a large fraction of the electrons in the flaring region become a 
part of the nonthermal spectrum, with a resulting energy content comparable to that of the 
magnetic field 5,6 . 

Mechanisms for particle acceleration have been explored and compared in a variety of 
papers, e.g.' -13 . Several authors 9,14,15 have examined acceleration by electric fields parallel 
to the local magnetic field (-Ej|)- However, parallel electric fields are typically localized near 
reconnection X-lines and separatrices, which limits the number of electrons that can be 
accelerated through this mechanism. 

Drake et al. 16 proposed a mechanism whereby charged particles gain energy as they reflect 
from the ends of contracting magnetic islands, a process analogous to the first-order Fermi 
acceleration of cosmic rays. This mechanism operates wherever there are contracting field 
lines in a reconnection region, and therefore develops during single X-line reconnection or 
as magnetic islands merge 12,13,16-18 . This mechanism is therefore volume filling and can 
accelerate a large number of particles. 

In a recent article 19 , we developed a method for calculating electron acceleration due to 
three fundamental mechanisms: parallel electric fields E n, betatron acceleration associated 
with conservation of the magnetic moment, and Fermi reflection due to the relaxation of 
curved magnetic field lines. We found that Fermi reflection dominated in reconnection where 
the magnetic fields are roughly antiparallel (see also Guo et al. 20 ), whereas in guide field 
reconnection Fermi reflection and E\\ were both important drivers of particle acceleration. 
However, we did not address the scaling of each mechanism with particle energy. This 
is important for determining the mechanism responsible for producing the most energetic 
particles. 

Studies of particle acceleration in reconnection have primarily been based on 2D simula¬ 
tions, in which accelerated particles are typically localized near the X-line, along magnetic 
separatrices and within magnetic islands 9,20 . There are some observations with small ambi- 
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ent guide fields 21-23 that support such a picture. A notable exception are Wind observations 
in which energetic electrons up to 300 keV are seen for more than an hour in an extended 
region around the reconnection region 1 . These observations correspond to reconnection with 
a strong guide held. 

Two-dimensional simulations impose limitations on the magnetic topology as well as the 
available spectrum of instabilities. In the presence of an ambient guide held, reconnection in 
3D can become turbulent as a result of the generation of magnetic islands along separatrices 
and adjacent surfaces 24,25 . While test particle trajectories in MHD helds have been used 
to explore acceleration in such systems, 26,27 the absence of feedback of energetic particles 
on the reconnection process in such models limits their applicability to real systems. Re¬ 
cent 3D studies of kinetic reconnection examined particle acceleration in electron-positron 
plasmas 20,28 . However, these studies focused on relativistic regimes where the magnetic en¬ 
ergy per particle exceeds the rest mass energy and included no ambient guide held. Hence, 
the impact of complex 3D magnetic helds on particle acceleration remains an open topic. 

Here, we explore magnetic reconnection in 3D systems with a strong guide held (i.e. 
one sufficiently strong to magnetize the electrons and to play an important role in pressure 
balance), which is the most generic form of reconnection in space and astrophysical plasmas. 
We find that the efficiency of particle acceleration is greatly increased compared to that in 
2D systems (where d/dz = 0, as in Dalilin et al. 19 ). We show that this occurs because 
the complex 3D magnetic helds enable the most energetic particles to continually access 
volume-filling acceleration sites rather than being confined to a single magnetic island that 
no longer accelerates particles once it has fully contracted. We also examine the energy 
dependence of the dominant E\\ and Fermi acceleration mechanisms, and fold that Fermi 
reflection is the primary accelerator of the energetic electrons. 

We explore particle acceleration via simulations using the massively parallel 3D particle- 
in-cell (PIC) code p3d 29 . Particle trajectories are calculated using the relativistic Newton- 
Lorentz equation, and the electromagnetic helds are advanced using Maxwell’s equations. 
The time and space coordinates are normalized, respectively, to the proton cyclotron time 
D" 1 = rriic/eB and inertial length di = c/co pi . The grid cell width is d e / 4, where d e = 
diy/m e /rrii is the electron inertial length. The time step is dt = O.OID^ 1 = 0.25D" 1 , where 
Vt ce = ( rrii/m e )Vt C i is the electron cyclotron frequency. 

We focus on a 3D simulation with dimensions L x x L y x L z = 51.2 di x 25.6 di x 25.6 di 
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and an analogous 2D simulation with L x x L y = 51.2 di x 25.6 d{. These simulations use an 
artificial proton-to-electron mass ratio mi/m e = 25 in order to reduce the computational 
expense. Simulations with differing mass-ratios and domains are presented to demonstrate 
the generality of the results. 

All simulations are initialized with a force-free configuration and use periodic boundary 
conditions. This is chosen as the most generic model for large-scale systems such as the 
solar corona where the density jump between the current layer and upstream plasma is 
not expected to be important. The magnetic field is given by: B x = B 0 tanh(y/it; 0 ) and 
B z = \]2 Bq — B 2 r . corresponding to an asymptotic guide field B z oo = B x oo = B 0 . We 
include two current sheets at y — L y / 4 and 3L y /4 to produce a periodic system, and 
w o = 1.25d e . This initial configuration is not a kinetic equilibrium, which would require a 
temperature anisotropy 32 , but is in pressure balance. 

The 3D simulations use at least 50 particles per cell for each species, and the 2D simula¬ 
tions use 1600 particles per cell. The initial electron and proton temperatures are isotropic, 
with T e = Ti = 0.25mjC^, and the initial density n 0 and pressure p are constant so that 
/3 = 8np/B 2 = 0.5. The speed of light is c = 3c ,4 y/rrii/m e , where ca = B 0 /^/4:irmin 0 . 

Reconnection develops from particle noise via the tearing instability, generating inter¬ 
acting flux ropes that grow and merge until they reach the system size at tfl ci ~ 50. The 
macroscopic evolution of the 2D and 3D systems is similar at this point, though the 2D 
simulation has released roughly 15% more magnetic energy. Fig. 1 shows an isosurface of 
one component of the electron current density J ez at tfl C i = 50 in the 3D simulation. The 
current exhibits filamentary structure that develops from instabilities with k z % 0 that are 
prohibited in 2D reconnection simulations 25 . 

In Fig. 2, energy spectra are shown for a variety of simulations in 2D and 3D with differing 
domain sizes and mass ratios. The spectra reveal significant electron acceleration in both 
2D and 3D simulations. However, the 3D simulations produce a much greater number of 
energetic particles: the fraction of electrons with energy exceeding 0.5m e c 2 is roughly an 
order of magnitude larger than in the 2D simulations (~ 2 x 10 " 4 vs. ~2x 1(T 5 ). The 
separation between the spectra in the 2D and 3D systems is greatest in the largest simulation 
domain (compare panels (a) and (c)) suggesting that this difference would be even greater 
for larger systems. Since the magnetic energy dissipation is typically slightly larger in the 2D 
systems, the increased energetic electron production in the 3D systems is due to enhanced 
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acceleration efficiency rather than an increase in the total energy imparted to the plasma. 
These results are insensitive to the mass ratio (panel (b)). The discussion that follows will 
focus on the results of the simulations shown in panel (c). 

The spatial distribution of the most energetic particles (shown in the left-hand panels of 
Fig. 3) also differs between the 2D and 3D simulations: these particles occupy narrow bands 
well inside the islands in the 2D simulation, but are distributed throughout the reconnecting 
region in the 3D simulation. In the 2D system, the reconnected held lines form closed loops 
(islands) that trap particles. The stochastic structure of the magnetic held in the 3D system, 
however, allows held-line-following particles to wander throughout the chaotic reconnecting 
region 33 . Surfaces of section from held line tracing reveal that the region of stochastic 
magnetic held in the 3D simulation roughly matches the distribution of energetic electrons. 
The distribution of the energetic particles in the 3D simulation is broadly consistent with 
Wind magnetotail observations where energetic electrons are seen for more than an hour 
and therefore must occupy a large region 1 . 

In order to examine the mechanisms responsible for accelerating these particles, we as¬ 
sume a guiding-center approximation relevant for a strong guide held 34 . In this limit, the 
evolution of the kinetic energy e of a single electron can be written as: 

J t = Q e \\ v \\ + ^ + u e • + 'ymevf K • k) (1) 

where E\\ = E-b is the parallel electric held, /i = m e ^ 2 vj_/2B is the magnetic moment, u B is 
the E x B velocity corresponding to the advection of the magnetic held, and k = b-Vb is the 
magnetic curvature. The velocity components parallel and perpendicular to the magnetic 
held are wy and v±, respectively; 7 is the relativistic Lorentz factor, and b is the unit vector 
in the direction of the local magnetic held. 

The hrst term on the right-hand-side of the equation corresponds to acceleration by 
parallel electric helds, which are typically localized near the reconnection X-line and along 
separatrices. The second term corresponds to betatron acceleration, which is typically negli¬ 
gible in reconnection 19 . The last term corresponds to reflection of particles from contracting 
magnetic held lines, a type of hrst-order Fermi acceleration 13,16,17 . This occurs where tension 
is released as magnetic helds advect in the direction of magnetic curvature (u B ■ n > 0). 
Both E || and Fermi reflection change the parallel energy of the particles, while betatron 
acceleration changes the perpendicular energy. 
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Equation (1) reveals that the acceleration mechanisms have different scalings with the 
particle energy: the Fermi reflection term is second-order in the parallel velocity, whereas 
the parallel electric field term is only first-order. Panel (d) of Fig. 2 shows the average 
acceleration per particle for both E\\ and Fermi reflection in the 3D simulation at tfl ci = 50. 
The bulk thermal electrons (low energies) are primarily accelerated by E n, whereas Fermi 
reflection is more important at high energies, consistent with the energy scaling of Eq. (1). 

The spatial distribution of the Fermi reflection term for the most energetic electrons 
(> 0.5m e c 2 ) is shown on the right-hand side of Fig. 3. While acceleration occurs throughout 
the reconnection exhaust in the 3D simulation, in 2D the acceleration is limited to narrow 
bands near the cores of magnetic islands. This contrast suggests that the stochastic 3D 
field structure allows the electrons to have greater access to the acceleration regions where 
magnetic energy is being released. 

To explore the reason for enhanced acceleration in the 3D system, we examine the trajec¬ 
tories of the 750 most energetic electrons in each simulation. A typical trajectory from the 
2D simulation is shown in the top left panel in Fig. 4. The electron begins in the tail of the 
electron distribution with kinetic energy e ~ 0.4m e c 2 . The electron streams along a field line 
outside the reconnection region before accelerating at an X-line near x ~ 50 and becoming 
trapped in an island. The electron bounces several times inside this island, accelerating up 
to e ~ 0.8m e c 2 . By this point, the field line the electron is following has released its tension, 
so acceleration ceases even as the electron continues to bounce. 

The top right panel of Fig. 4 shows a typical electron trajectory from the 3D simulation. 
The electron moves throughout the reconnecting domain, and does not become trapped on 
stagnant field lines in island cores as in the 2D simulation. It instead gains energy at many 
different acceleration sites, eventually reaching a maximum energy of e ~ 1.15m e c 2 . The 
Supplemental Material 38 includes videos of both particle trajectories in the time-evolving 
field of the simulation. 

The electron trajectories shown here are generic for their respective simulations. Though 
the acceleration details differ, all of the electrons in the 2D simulation are confined to 
single islands, whereas no electrons in the 3D simulation show significant trapping. The 
bottom panels of Fig. 4 show the distribution of the particle displacement | Ax| = \x(Q C it = 
50) — x{VL ci t = 25)| for the 750 most energetic particles in each simulation (the choice of 
Q ci t = 25 as the earliest time eliminates free streaming along unreconnected field lines 
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before islands develop). The average displacement of the energetic electrons in the 3D 
system is nearly an order of magnitude greater than that in the 2D simulation, underscoring 
a fundamental difference in the particle trajectories of the two systems. 

It has been shown previously that the development of pressure anisotropy with Py ^ P± 
causes the cores of magnetic islands to approach firehose marginal stability, where the tension 
driving magnetic reconnection ceases, thereby throttling reconnection. A full treatment of 
the feedback from particle acceleration (e.g. Drake et al. 18 ) is outside the scope of this paper. 
However, we do find that large anisotropies persist in the 3D simulations (Supplemental 
Material 38 ), so the turbulent dynamics do not appear to significantly isotropize the pressure. 
The firehose instability is not triggered in these simulations, likely due to the strong ambient 
guide field - there is insufficient magnetic free energy to reach /3 ~ 1. In cases with weaker 
guide fields the firehose feedback is likely to be more important 35 . 

The electron spectra in both simulations do not assume a power law as is frequently 
observed in nature. This is due in part to the limited energy gain possible in the modest¬ 
sized 3D simulation presented here. Previous 2D simulations have shown the total energy 
gain is greater in larger systems 19 . An additional issue is that these simulations have periodic 
boundary conditions so no particles are lost from the system. Solar observations suggest 
that electrons are confined in regions of energy release in the corona 5 . The mechanism for 
confinement remains an open issue, although both magnetic mirroring and double layers are 
possible mechanisms 36 . On the other hand, it has been suggested that the development of 
a power law requires a loss mechanism in addition to an energy drive 18 . However, recent 
electron-positron simulations 20,28 suggest that power-law spectra may still develop in the 
absence of a loss mechanism. The set of conditions under which power law spectra form in 
kinetic reconnection simulations remains an open issue. 

A limitation of the present simulations is the use of an artificial mass ratio, which reduces 
the separation between proton and electron scales. In order for an electron to access multiple 
acceleration sites as we observe in our simulations, its characteristic velocity must exceed that 
of the macroscopic flows associated with reconnection, which is controlled by the protons. 
The results of a simulation performed with mj/m e = 100 were not qualitatively different 
from the m,i/m e = 25 simulation, suggesting we have achieved a sufficient separation of scales 
(see Fig. 2b). The absence of such a separation of scales may explain why this behavior is 
not observed in electron-positron simulations 20 . 
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FIG. 1. Isosurface of J ez at tQ c i = 50. The isosurface level is 60% of the maximum current density 
(a 2D slice of the same quantity is shown on the bottom). 

The plasma parameters used in these simulations are not tuned to a specific astrophysical 
system. For example, the ratio of the speed of light to the electron Alfven velocity (c/ CA e = 3) 
is smaller than is typical for either the magnetotail 4 or the solar wind 3 '. However, so long as 
c/cAe > 1, the important reconnection dynamics are not relativistic. The initial temperature 
(T e0 = rriic\/8) is much larger than in the either the corona 5 or the magnetotail 4 . As long as 
the electron thermal velocity Vt e o ca, electrons circulate in islands faster than the islands 
contract and a preheating mechanism is not needed for Fermi reflection to function. 
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FIG. 2. (a-c) Global electron energy spectra. The 3D simulation dimensions L x x L y x L z are: 
(a) 102.4 x 51.2 x 25.6, (b) 51.2 x 25.6 x 12.8, (c) 51.2 x 25.6 x 25.6. Dotted lines indicate initial 
spectra, solid lines in (b-d) correspond to t = 50. Dashed and solid lines in (a) correspond to t = 50 
and t = 125, respectively, (d) Average electron energization rate vs. energy for the 3D simulation 
shown in (c). 


Simulations were carried out at the National Energy Research Scientific Computing Center. 
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